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Abstract 

Available  humidity  sensing  techniques  are  often  intrusive,  and  of  limited  practical  interest  for  real-time  control  applications  due  to  their 
cost,  size,  and  inadequate  response  time  and  accuracy.  In  this  study,  we  present  a  novel  method  for  estimation  of  PEM  fuel  cell  humidity  by 
exploiting  its  effect  on  cell  resistive  voltage  drop.  This  voltage  loss  is  discerned  from  mass  transport,  concentration,  activation  losses  and  open 
circuit  voltage  by  a  well-known  fuel  cell  voltage  model.  The  proposed  scheme  makes  use  of  measurements  of  voltage,  current,  temperature, 
and  total  pressure  values  in  the  anode  and  cathode.  It  also  incorporates  dynamic  estimators  for  hydrogen  and  oxygen  partial  pressures,  adapted 
from  [M.  Arcak,  H.  Gorgun,  L.M.  Pedersen,  S.  Varigonda,  A  nonlinear  observer  design  for  fuel  cell  hydrogen  estimation,  IEEE  Trans.  Control 
Syst.  Technol.  12  (1)  (2004)  101-110].  The  membrane  resistance  thus  obtained  is  then  used  to  estimate  membrane  water  content  following 
functional  characterizations  presented  in  [T.E.  Springer,  T.A.  Zawodzinski,  S.  Gottesfeld,  Polymer  electrolyte  fuel  cell  model,  J.  Electrochem. 
Soc.  138  (8)  (1991)  2334-2342].  Experiments  with  this  estimation  technique,  performed  at  the  Connecticut  Global  Fuel  Cell  Center,  are 
presented  and  discussed. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Polymer  electrolyte  membrane  (PEM)  fuel  cells  are  envi¬ 
sioned  to  be  the  future  choice  for  portable  power,  trans¬ 
portation,  and  combined  heat  and  power  systems.  They  offer 
advantages  to  other  fuel  cell  types,  including  relative  sim¬ 
plicity  of  their  design,  and  their  ability  to  operate  at  low 
temperatures.  In  PEM  fuel  cells  the  membrane  must  be  suf¬ 
ficiently  hydrated  because  its  conductivity  depends  critically 
on  the  humidity  level.  Too  little  water  causes  membrane 
drying,  which  increases  the  ionic  resistance,  and  exacer¬ 
bates  the  voltage  drop  due  to  ohmic  losses.  Too  much  water 
causes  “flooding”,  that  is  blocking  of  porous  passages,  which 
reduces  the  transport  rate  of  reactants  to  the  catalyst  site. 
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An  obstacle  to  active  control  of  membrane  water  content 
is  the  lack  of  adequate  tools  for  monitoring  humidity  in  the 
fuel  cell.  The  cost  and  size  of  existing  humidity  sensors  are 
prohibitive  for  in  situ  measurements.  Their  accuracy  is  also 
impaired  for  high  relative  humidity  levels  [3] ,  at  which  a  PEM 
fuel  cell  must  operate.  Other  measurement  techniques,  such 
as  [4]  which  employs  gas  chromatography,  require  extractive 
sampling  and,  thus,  are  slow  and  intrusive. 

In  this  paper,  we  present  an  estimation  scheme  for  the 
membrane  water  content.  Our  idea  is  to  first  estimate  the 
membrane  resistance  and,  next,  to  use  available  characteri¬ 
zations  of  membrane  resistance  as  a  function  of  water  content, 
such  as  those  in  Ref.  [2]  for  Nation  112  membranes.  To 
accomplish  the  first  step  of  membrane  resistance  estima¬ 
tion  we  calculate  the  ohmic  voltage  loss  from  voltage  and 
current  measurements,  and  from  the  well-known  fuel  cell 
voltage  model  [5]  which  accounts  for  other  voltage  loss  terms. 
Because  the  open-circuit  voltage  component  of  this  model 
depends  on  partial  pressures  of  hydrogen  and  oxygen,  which 
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are  unavailable  for  measurement,  we  estimate  them  from  a 
variant  of  the  hydrogen  and  oxygen  observers  developed  in 
Ref.  [1].  We  wish  to  emphasize  that  our  estimation  scheme  is 
applicable  when  the  current  is  nonzero,  because  it  relies  on 
the  resistive  voltage  drop  which  is  induced  by  the  current. 

A  different  approach  to  humidity  estimation  is  presented  in 
Ref.  [6] ,  where  the  authors  employ  open-loop  observers  based 
on  lumped  dynamic  models  for  anode  and  cathode  relative 
humidity  values.  This  approach  differs  from  our  voltage- 
based  estimation  because,  first,  it  does  not  make  use  of  the 
voltage  output  of  the  humidity  model  and,  second,  it  assumes 
open-circuit  conditions.  In  contrast,  we  rely  on  the  voltage 
output  and  its  static  relation  to  water  content,  and  do  not 
employ  a  dynamic  humidity  model. 

The  paper  is  organized  as  follows:  In  Section  2  we  review 
the  fuel  cell  voltage  model  and  the  characterization  of  mem¬ 
brane  resistance  as  a  function  of  its  water  content.  The 
estimation  algorithm  is  detailed  in  Section  3,  followed  by 
experimental  results  and  their  interpretations  in  Section  4. 
Remaining  research  tasks  and,  in  particular,  a  discussion  of 
how  the  estimation  design  can  be  modified  to  account  for 
flooding  conditions,  are  presented  in  Section  5. 


2.  Overview  of  the  voltage  model 


The  theoretical  (reversible)  voltage  of  a  fuel  cell  is  given 
by  (see  e.g.  ([5],  Chapter  2)): 


E 


A  g°{  RT 
IF  +  IF 


1 

H —  In 
2 


where  Agj?  <0  is  the  change  in  molar  Gibbs  free  energy 
of  formation  at  standard  pressure,  Pq,  R  the  universal  gas 
constant,  F  the  Faraday  constant,  T  the  cell  temperature, 
ln(-)  denotes  the  natural  logarithm,  and  and  repre¬ 
sent  average  values  of  hydrogen  and  oxygen  partial  pressures 
across  the  anode  and  cathode  channels,  respectively.  As  dis¬ 
cussed  in  Ref.  [1]  (Section  III),  a  good  approximation  for 
these  average  values  is  the  arithmetic  mean  of  the  inlet  and 
exit  partial  pressures.  However,  in  some  situations  the  electro- 
osmotic  drag  may  not  be  constant  along  the  channel,  thus 
affecting  the  accuracy  of  this  approximation. 

The  operational  voltage  of  the  fuel  cell  differs  from  its 
theoretical  value,  E\  due  to  the  activation  voltage  drop ,  Va; 
the  ohmic  voltage  drop  V0hm;  the  mass  transport  loss  Vmass, 
as  detailed  in  Ref.  [5]  (Chapter  3).  For  a  stack  of  n  cells  for 
which  the  losses  are  identical,  the  net  voltage  is: 


Vst  —  n(E  Va  Vohm  Vmass)- 


The  cumulative  effect  of  the  voltage  loss  terms  in  (2)  is 
visible  from  the  polarization  curve  in  Fig.  1,  which  gives 
a  plot  of  the  cell  voltage  as  a  function  of  the  current  den¬ 
sity ;  that  is,  total  current  I  divided  by  effective  membrane 


Fig.  1 .  Voltage  of  a  single  cell  (V)  vs.  current  density,  i  (A  cm-2).  The  voltage 
model  of  Section  2  (continuous  curve)  matches  the  discrete  data  points. 
Running  conditions  are:  Tce\\  =  60  °C,  H2  flow  =  1  slpm,  air  flow  =  1.5  slpm, 
humidification  temperature  =  60  °C. 


area,  Am: 

I 


Other  voltage  losses  due  to  fuel  crossover  and  internal 
currents  ([5],  Section  3.5)  are  not  included  in  (2)  because 
their  effect  is  considerable  only  at  very  low  current  densities. 

The  activation  voltage  drop  is  the  voltage  lost  in  driving 
the  chemical  reactions  on  the  surface  of  the  electrodes.  It  is 
given  by  the  well-known  empirical  formula  ([5],  Section  3.4): 


Fa  = 


RT 
2 aF 


i  >  i0, 


where  a  is  the  charge  transfer  coefficient,  and  io  is  the 
exchange  current  density,  which  depends  on  the  temperature, 
pressure,  the  type  of  catalyst  and  its  specific  surface  area  and 
loading  [7] .  Likewise,  the  mass  transfer  voltage  drop  is  given 
by  ([5],  Section  3.7): 


Vmass  —  B  In  M  .  ]  i  <  fiim?  (5) 

\  fiim  / 

where  B  is  a  constant  that  depends  on  the  fuel  cell,  and  its 
operating  conditions  and  fiim  is  the  limiting  current  density , 
estimated  from  experimental  data. 

Finally,  the  ohmic  voltage  drop  is 


Fohm  —  (  R\\\  +  T )/, 


where  Rm  is  the  membrane  resistance,  and  r  represents  the 
sum  of  other  components,  such  as  resistance  through  electri¬ 
cally  conductive  components  of  the  fuel  cell,  including  con¬ 
tact  resistance,  and  resistance  from  the  electrodes.  Because 
measurement  techniques  are  available  to  distinguish  between 
the  components  of  total  resistance  [8],  for  our  design  in  the 
next  section  we  assume  that  r  is  known. 

The  membrane  resistance,  Rm ,  depends  critically  on  the 
membrane  water  content ,  which  is  defined  as  the  ratio  of  the 
number  of  water  molecules  to  the  number  of  charge  sites. 
Indeed,  as  shown  in  Ref.  [2],  water  content  affects  the  mem- 
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brane  resistance  via 


Am(0.00514Am  -0.00326) 


x  exp 


f  1  1 

1268  - 


T  303 


where  tm  and  Am  are  the  membrane  thickness,  and  area, 
respectively.  The  study  in  Ref.  [2]  also  relates  the  membrane 
water  content  to  the  relative  humidity ,  0m  (the  ratio  of  water 
partial  pressure  to  saturated  vapor  pressure)  that  the  mem¬ 
brane  is  exposed  to  by  the  empirical  formula: 

V  =  0.043  +  17.81</>m  -  39.85</4  +  36.0^.  (8) 


3.  Estimation  of  water  content 


To  estimate  the  water  content,  Am,  our  approach  is  to  first 
estimate  Rm  and,  next,  to  invert  the  function  (7)  to  obtain 

/V 

the  estimate  km.  For  the  estimation  of  Rm  we  rely  on  (6),  in 
which  /  is  available  for  measurement  and  r  is  known.  Because 
Vohm  is  not  measured  separately,  it  is  to  be  calculated  from 
the  measurement  of  the  net  voltage  Vst,  by  substituting  in  (2) 
the  values  of  E ,  Va,  and  Vmass,  calculated  from  (1),  (4),  and 
(5),  respectively. 

This  algorithm  should  not  be  applied  for  zero  or  small  val¬ 
ues  of  the  current  /,  because  the  computation  of  Rm  from  (6) 
involves  division  by  I.  It  is  reliable  when  the  fuel  cell  oper¬ 
ates  in  the  linear  region  of  the  polarization  curve  because, 
then,  the  effect  of  humidity  on  the  open-circuit  voltage,  not 
modeled  in  (1),  and  losses  due  to  fuel  crossover  and  internal 
currents,  not  accounted  for  in  (2),  are  indeed  negligible  com¬ 
pared  to  the  ohmic  voltage  drop.  Since  the  linear  region  of  the 
polarization  curve  is  the  desired  regime  of  fuel  cell  operation, 
this  restriction  does  not  impair  the  practical  relevance  of  our 
algorithm. 

The  main  difficulty  in  this  algorithm  is  the  calculation 
of  the  open-circuit  voltage,  E ,  from  (1),  which  relies  on  the 
unmeasured  hydrogen  and  oxygen  partial  pressures,  and 
Pq?  .  To  overcome  this  difficulty  we  employ  a  variant  of  the 
hydrogen  and  oxygen  observer  developed  in  Ref.  [1].  To  this 
end  we  denote  by  pu2  the  exit  partial  pressure  of  hydrogen 
in  the  anode  channel  and  obtain,  from  the  Ideal  Gas  Law,  the 
lumped  dynamic  model  [9] : 


RT 

Vol^ 


1  PH2in  F  . _ 1  PH2  f 

Pain  “  Mao  Pa  ! 


where  /?h2  denotes  the  time  derivative  of  pu2 •  Vola  is  the 
anode  volume,  Fa  the  total  anode  pressure,  pn2m  and  Fain, 
respectively,  the  hydrogen  partial  pressure  and  the  total  pres¬ 
sure  at  the  inlet  of  the  anode,  Fai  and  Fao  the  total  inlet  and 
outlet  mass  flows,  and  Ma i  and  Mao  represent  average  molec¬ 
ular  weights  of  the  inlet  and  exit  gas  streams.  The  bracketed 
expression  in  (9)  calculates  the  molar  rate  of  change  of  hydro¬ 
gen.  In  particular,  the  first  term  represents  the  increase  in 
hydrogen  mole  number  due  to  the  inlet  flow,  the  second  term 


is  the  decrease  due  to  the  exit  flow,  and  the  third  term  is  the 
consumption  due  to  the  fuel  cell  reaction. 

For  observer  design  we  assume  that  the  total  pressures  Fain 
and  Fa,  and  the  mass  flows  Fa i  and  Fao,  are  known.  The  mass 
flows  can  either  be  measured  with  flow  meters,  or  estimated 
from  pressure  differences  via  orifice  equations,  such  as  those 
used  in  Ref.  [1].  The  variables  Ma i  and  pnrm  which  appear  in 
the  coefficient  of  Fa i  in  (9)  are  also  assumed  to  be  available. 
This  assumption  is  meaningful  when  the  gas  composition 
of  the  inlet  stream  is  known,  or  estimated  with  a  separate 
observer  for  the  fuel  reformer  as  in  Ref.  [10].  Our  observer 
mimics  Eq.  (9)  and  produces  an  estimate,  pn2 ,  for  the  exit 
partial  pressure  pn2  from: 


For  the  unknown  Mao,  we  employ  the  value 


Pc i  Fh2 

P^~ 


(11) 


in  which  Mh2  is  the  molecular  weight  of  hydrogen,  and  8a  is 
an  average  value  for  the  molecular  weights  of  other  gases  at 

/V 

the  anode  exit.  If  pure  hydrogen  is  supplied  to  the  anode,  Mao 
must  be  taken  to  be  Mh2  .  Unlike  the  design  in  Ref.  [1]  which 
employed  a  constant  average  value  for  Mao,  the  new  observer 
(10)  and  (11)  is  nonlinear  in  pn2  due  to  the  dependence  of 

/V 

Mao  on  pu2 .  However,  because  hydrogen  is  lighter  than  the 
other  molecules;  that  is,  Mh2  <  <$a  in  (11),  the  right-hand 
side  of  (10)  is  a  monotone  decreasing  function  of  pu2 ,  which 
guarantees  observer  convergence  as  proven  in  Ref.  [11]. 

Using  a  similar  model  for  the  cathode,  and  denoting  by 
po2  the  estimate  of  po2  we  obtain  the  oxygen  observer 


(12) 


(13) 


where  Volc  is  the  cathode  volume,  Pc  the  total  cathode  pres¬ 
sure,  po2in  and  Fcin  the  oxygen  partial  pressure  and  the  total 
pressure  at  the  inlet  of  the  cathode,  Fc{  and  Fco  the  total 
inlet  and  outlet  mass  flows,  Mc{  and  Mco  represent  average 
molecular  weights  of  the  inlet  and  exit  gas  streams,  Mo2  the 
molecular  weight  of  oxygen  and  8C  is  an  average  value  for 
the  molecular  weights  of  other  gases  in  the  exit  stream. 

To  estimate  the  membrane  water  content  km,  we  first  esti¬ 
mate  the  theoretical  voltage  E  by  obtaining  pn2  and  po2  from 
the  observers  (10),  (1 1)  and  (12),  (13),  and  by  substituting  in 
(1)  the  average  values: 


„av  _  Ph2  +  Pli2m  „av  _  P02  +  P02in 

p n2  -  5  Po2  -  9 


(14) 


Next,  we  calculate  V0hm  from  (2),  using  our  estimate  of 
F,  as  well  as  the  measurement  of  voltage  Vst,  and  the  calcu¬ 
lation  of  Va  and  Vmass  from  (4)  and  (5).  We  then  obtain  the 
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/v 

membrane  resistance  estimate  Rm  from  (6)  by  using  the  cur¬ 
rent  measurement  /,  and  by  subtracting  r.  Finally,  we  use  this 

✓V  /V  /V 

Rm  to  estimate  the  water  content  from  (7),  and  obtain  </>m 
from  the  unique  real  root  of  the  polynomial  (8),  given  by: 

4>m  =  0.369  +  (S  +  VQ3  +  S2)l/ 3  +  (S  -  VQ3  +  S2)1 3 

(15) 

where  S  =  —0.0416  +  Am/72  and  Q- 0.0288. 

4.  Experimental  results 

We  now  present  experimental  results  for  the  algorithm 
developed  in  the  previous  section.  These  experiments  have 
been  performed  at  the  Connecticut  Global  Fuel  Cell  Center, 
with  Fuel  Cell  Technologies  Inc.  hardware  and  a  3  M  mem¬ 
brane  electrode  assembly.  The  active  area  is  Am  =  50  cm2,  and 
the  membrane  thickness  is  tm  =  0.0051  cm.  Both  anode  and 
cathode  flow  channels  are  four  path  serpentine.  Hydrogen 
and  air  are  supplied  from  hydrogen  and  air  cylinders.  Hydro¬ 
gen  and  air  flows  are  regulated  and  humidified  by  a  Teledyne 
Medusa  RD  fuel  cell  test  station.  A  Scribner  Associates  890C 
electronic  load  bank  is  used  to  generate  the  load  current  pro¬ 
file.  The  cell  resistance  is  measured  by  the  current  interrupt 
technique  which  is  integrated  to  the  load  bank  and  controlled 
by  software. 

We  first  compare  the  voltage  model  of  Section  2  to  the 
experimental  polarization  curve.  To  obtain  this  curve  we  set 
the  hydrogen  flow  rate  to  1  standard  liter  per  minute  (slpm), 
and  the  air  flow  rate  to  1 .5  slpm.  The  hydrogen  and  air  humid¬ 
ification  temperatures  were  set  to  60  °C.  To  prevent  conden¬ 
sation,  we  set  the  line-heaters  for  hydrogen  and  oxygen  to 
70  °C.  We  observed  the  open  circuit  voltage  for  each  cell,  then 
gradually  increased  the  current  up  to  40  A  (800  mA  cm-2  at 
H2  stoichiometry  around  3),  and  reduced  it  back  to  zero.  The 
resulting  curve  is  given  in  Fig.  1,  where  the  discrete  data 
points  are  obtained  from  step  changes  in  the  current  every 
3  min.  We  obtained  a  close  match  to  these  data  points  using 
the  voltage  model  of  Section  2  (continuous  curve  in  Fig.  1) 
with  B  =  RT/2F  and  iyim  =  1 .4  A  cm-2  in  (5),  a  =  0.5  in  (4),  and 
with  z'o  obtained  from  a  temperature-  and  pressure-dependent 
formula  in  Ref.  [7],  in  which  we  used  the  reference  value 
*Qef  =  3.4  x  10-11  A  cm-2,  resulting  in  apparent  exchange 
current  density  of  z'o  =  8.8  x  10-8  A  cm-2. 

To  test  our  humidity  estimation  algorithm,  we  induced 
membrane  drying  by  increasing  the  cell  temperature  above 
the  humidification  temperature.  Specifically,  the  cell  tem¬ 
perature  was  raised  from  60  to  80  °C  while  the  humidifier 
temperature  was  kept  at  60  °C.  This  causes  the  air  to  absorb 
water  from  the  membrane  as  it  passes  through  the  warmer 
cell.  The  resulting  increase  in  the  resistance  and  the  decrease 
in  the  cell  voltage  are  shown  in  Fig.  2  for  a  current  of  40  A. 
Due  to  the  absence  of  humidity  measurements  in  our  set-up, 
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Fig.  2.  Voltage  decrease  and  resistance  (x  10  2  £2)  increase  in  response  to 
drying,  which  starts  at  t=  600  s.  Solid  line  is  the  measured  resistance  and  the 

v\ 

dashed  line  is  the  sum  of  our  estimate  Rm,  and  r  =  0.344  mQ. 

resistance  to  that  obtained  from  the  current  interrupt  mea¬ 
surement  in  Fig.  2.  The  observer  calculations  are  presently 
carried  out  off-line,  from  filtered  experimental  data. 

At  fully  humidified  conditions  the  cell  performance  was 
in  steady  state  and  the  resistance  measurement  indicated 
1.33  m^  (65m£2cm2),  as  shown  in  Fig.  2.  The  model  (7) 
predicts  the  membrane  resistance  to  be  Rm  =  0.986  mQ 
(49m£2cm2),  which  means  r-  1.33-0.986  =  0.344  mQ 
(17m£2cm2).  lOm^cm2  of  this  may  be  attributed  to  the 
resistance  from  the  electrodes  [12].  Therefore,  the  electronic 
resistance  is  7  m£2  cm2,  which  includes  the  resistance 
through  the  electrically  conductive  components  of  the  fuel 
cell.  In  Ref.  [13],  the  electronic  resistance  is  reported  to  be 
up  to  three  times  larger  than  the  ionic  resistance  (membrane 
and  electrodes).  This,  of  course,  depends  on  the  cell  design. 

Under  drying,  the  cell  resistance  increased  at  a  rate  of 
0.0018  m^s-1.  The  electronic  resistance  should  be  fairly 
independent  of  the  drying  conditions,  i.e.,  membrane  water 
content.  Although  drying  of  the  membrane  results  in  reduced 
thickness,  in  a  well  designed  fuel  cell  stack,  the  changes  in 
the  membrane  thickness  should  be  compensated  by  a  stack 
compression  mechanism.  It  is  therefore  safe  to  assume  that 
only  the  ionic  resistance  changes  with  humidity  in  the  cell. 

The  resulting  observer  estimate  of  Rm  derived  from  the 
observance  of  the  cell  stack  potential,  upon  addition  of 
r-  0.344  m^  is  shown  with  a  dotted  curve  in  Fig.  2.  It 
achieves  reasonable  accuracy  with  respect  to  the  measured 
resistance  (solid  curve).  After  t=  1000  s,  however,  it  over¬ 
estimates  the  measured  values.  This  is  due  to  other  factors 
affecting  the  cell  potential  under  drying  conditions,  such  as 
the  loss  of  active  sites.  Indeed,  it  is  reported  in  Ref.  [14] 
that  drying  of  the  anode  catalyst  layer  not  only  increased  the 
membrane  proton-conduction  resistance  but  also  increased 
the  activation  overpotential  for  the  hydrogen  oxidation  reac¬ 
tion  due  to  a  decrease  in  the  number  of  active  sites  on  the 
anode  side.  A  similar  argument  may  be  applied  to  the  cath¬ 
ode  drying  scenario,  where  the  loss  of  active  sites  results  in 
a  change  of  apparent  exchange  current  density.  In  our  exper¬ 
iment,  the  apparent  exchange  current  density  changed  from 
/o  —  8.8  x  1CD8  A  cm-2  at  60  °C  with  fully  humidified  reac¬ 
tant  gases  to  io  =  1 .01  x  10-8  A  cm-2  at  80  °C  with  dry  gases, 
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Fig.  3.  Impedance  spectroscopy  of  the  cell  under  humidified,  and  dry  con¬ 
ditions. 

which  corresponds  to  about  60  times  lower  active  catalyst 
surface  area. 

This  explanation  is  corroborated  by  impedance  spec¬ 
troscopy  results  in  Fig.  3  where  the  left  lobe  corresponds 
to  humid  conditions  (100%  relative  humidity)  while  the 
right  lobe  corresponds  to  dry  conditions  (40%).  The  mea¬ 
surements  were  collected  over  a  frequency  range  of  0. 1  Hz 
to  10  kHz  with  amplitude  of  10  mV  using  a  Solartron  Inc. 
electrochemical  interface  and  frequency  response  analyzer. 
The  left  intersection  with  the  real  axis  is  the  high  frequency 
measurement,  which  represents  membrane,  contact  and  elec¬ 
tronic  resistances,  while  the  right  intersection,  low  frequency 
measurement,  includes  charge  transfer  resistance  in  addi¬ 
tion.  Under  drying  conditions  the  left  intersection  moves  to 
the  right  indicating  an  increase  in  membrane  resistance.  The 
width  of  the  lobe,  however,  also  increases  which  means  that 
the  charge  transfer  resistance  has  also  increased. 

We  wish  to  emphasize,  however,  that  the  main  interest  in 
our  estimation  algorithm  is  for  use  in  a  feedback  controller 
that  regulates  the  humidity  around  100%  and,  thus,  the  accu¬ 
racy  achieved  within  the  first  1 100  s  in  Fig.  2  is  satisfactory 
for  appropriate  control  action  to  be  taken. 

/V 

Using  the  estimated  Rm  in  Fig.  2,  we  obtain  from  Eqs.  (7) 

/V  /\ 

and  (15)  the  estimates,  Am  and  <pm,  as  in  Fig.  4.  It  is  clear 


Fig.  4.  Water  content  estimate  (solid  line),  and  relative  humidity  estimate 

A 

0m  (dashed  line),  scaled  by  10,  obtained  via  Eqs.  (7)  and  (15)  from  the 
estimate  in  Fig.  2. 


/v 

that  the  membrane  was  intensively  drying  with  Am  dropping 
below  4.  Relative  humidity  of  the  gases,  once  they  reach  the 
stack  operating  temperature,  was  changing  from  100%  at  cell 
temperature  of  60  °C,  to  64%  at  70  °C,  to  42%  at  80  °C.  Cor¬ 
responding  relative  humidities  at  the  stack  outlet  (both  anode 
and  cathode)  were  100%  at  both  60  and  70  °C,  while  at  80  °C 
it  was  69%.  The  average  relative  humidity  of  the  gases  at 
80  °C  was  about  50%,  which  is  also  confirmed  in  Fig.  4. 

Note  that  the  experiments  above  do  not  include  flooding 
conditions  which  would  also  cause  a  voltage  drop  at  constant 
current.  To  distinguish  between  flooding  and  drying,  a  cath¬ 
ode  pressure  drop  may  be  monitored  [15]  and  incorporated 
in  our  algorithm.  Pressure  drop  on  the  cathode  side  increases 
with  cell  flooding,  while  it  remains  unchanged  with  cell  dry¬ 
ing,  thus  clearly  distinguishing  between  the  two  phenomena 
[15-17].  Further  details  on  the  use  of  cathode  pressure  drop 
as  a  diagnostic  signal  may  be  found  in  Ref.  [15],  and  is  not 
repeated  here. 

5.  Conclusions 

Water  management  is  an  important  problem  for  PEM 
fuel  cell  operation.  In  this  paper  a  method  for  estimation 
of  membrane  water  content  has  been  developed  by  exploit¬ 
ing  its  effect  on  cell  resistive  voltage  drop.  This  method 
relies  on  measurements  of  voltage,  current,  temperature,  and 
total  pressure  in  the  anode  and  cathode.  It  also  incorporates 
observers  for  hydrogen  and  oxygen  partial  pressures  adapted 
from  Ref.  [1]. 

A  voltage  drop  at  constant  current  may  also  be  caused 
by  flooding — an  effect  completely  opposite  from  drying.  In 
order  to  avoid  an  incorrect  attribution  of  this  voltage  drop  to 
drying,  a  cathode  pressure  drop  may  be  monitored  [15]  and 
incorporated  in  our  algorithm.  Pressure  drop  on  the  cathode 
side  increases  with  cell  flooding,  while  it  remains  unchanged 
with  cell  drying,  thus  clearly  distinguishing  between  the  two 
phenomena  [15-17]. 

Dry  conditions  at  the  cell  inlet  and  outlet  mean  that  all 
of  the  water  generated  in  the  electrochemical  reaction  evap¬ 
orated,  unlike  in  steady  state  saturated  conditions  at  cathode 
inlet  and  outlet  when  all  of  the  product  water  left  the  stack  as 
liquid.  The  Gibbs  free  energy,  and  therefore  the  theoretical 
cell  potential,  are  lower  for  gaseous  product  water.  Addi¬ 
tional  effects  of  water  evaporation  on  local  temperature,  and 
therefore  on  saturation  conditions,  may  only  be  determined 
by  detailed  3D  modeling,  and  then  integrated  over  the  entire 
cell. 

An  additional  effect  of  cell  drying  on  activation  over¬ 
potential  has  been  detected  in  experiments.  This  effect  can 
be  compensated  for  in  our  algorithm  upon  further  modeling 
efforts  for  the  increase  of  activation  overpotential  relative  to 
the  increase  in  purely  resistive  voltage  drop.  Compensation 
for  the  additional  voltage  losses  brought  on  by  time  (MEA 
degradation,  poisoning,  etc.)  would  be  possible  with  further 
diagnostic  tools  that  indirectly  detect  such  losses. 


394 


H.  Gorgun  et  al.  /  Journal  of  Power  Sources  157  (2006)  389-394 


Acknowledgements 

The  authors  would  like  to  thank  Dr.  Alevtina  Smirnova 
of  the  Connecticut  Global  Fuel  Cell  Center  for  her  help  in 
interpreting  experimental  results.  The  work  of  the  second 
author  was  supported  by  NSF  under  grant  ECS-0238268. 

References 

[1]  M.  Arcak,  H.  Gorgun,  L.M.  Pedersen,  S.  Varigonda,  A  nonlinear 
observer  design  for  fuel  cell  hydrogen  estimation,  IEEE  Trans.  Con¬ 
trol  Syst.  Technol.  12  (1)  (2004)  101-110. 

[2]  T.E.  Springer,  T.A.  Zawodzinski,  S.  Gottesfeld,  Polymer  electrolyte 
fuel  cell  model,  J.  Electrochem.  Soc.  138  (8)  (1991)  2334-2342. 

[3]  R.S.  Glass  (Ed.),  Sensor  Needs  and  Requirements  for  Proton- 
Exchange  Membrane  Fuel  Cell  Systems  and  Direct-Injection 
Engines,  Lawrence  Livermore  National  Laboratory,  Applied  Energy 
Technologies  Program,  Livermore,  California,  2000. 

[4]  M.M.  Mench,  Q.L.  Dong,  C.Y.  Wang,  In  situ  water  distribution  mea¬ 
surements  in  a  polymer  electrolyte  fuel  cell,  J.  Power  Sources  124 
(2003)  90-98. 

[5]  J.  Larminie,  A.  Dicks,  Fuel  Cell  Systems  Explained,  second  ed., 
John  Wiley  &  Sons  Inc.,  Chicester,  2003. 

[6]  D.  McKay,  A.  Stefanopoulou,  Parametrization  and  validation  of  a 
lumped  parameter  diffusion  model  for  fuel  cell  stack  membrane 
humidity  estimation,  in:  Proceedings  of  the  2004  American  Control 
Conference,  Boston,  MA,  2004,  pp.  816-821. 

[7]  H.A.  Gasteiger,  W.  Gu,  R.  Makharia,  M.F.  Mathias,  Tutorial:  catalyst 
utilization  and  mass  transfer  limitations  in  the  polymer  electrolyte 
fuel  cell,  in:  The  2003  Electrochemical  Society  Meeting,  Orlando, 
FL,  2003. 


[8]  M.  Smith,  D.  Johnson,  L.  Scribner,  Electrical  test  methods  for  eval¬ 
uating  fuel  cell  ME  A  resistance,  Fuel  Cell  Mag.  (February /March) 
(2004)  15-17. 

[9]  J.  Pukrushpan,  A.G.  Stefanopoulou,  H.  Peng,  Control  of  Fuel  Cell 
Power  Systems:  Principles,  Modeling,  Analysis  and  Feedback  Design 
(Advances  in  Industrial  Control),  Springer,  2004. 

[10]  H.  Gorgun,  M.  Arcak,  S.  Varigonda,  S.  Bortoff,  Nonlinear  observer 
design  for  fuel  processing  reactors  in  fuel  cell  power  systems,  in: 
Proceedings  of  the  2004  American  Control  Conference,  Boston,  MA, 
2004,  pp.  845-849. 

[11]  M.  Arcak,  P.  Kokotovic,  Nonlinear  observers:  a  circle  crite¬ 
rion  design  and  robustness  analysis,  Automatica  37  (12)  (2001) 
1923-1930. 

[12]  F.N.  Buchi,  G.G.  Scherer,  Investigation  of  the  transversal  water 
profile  in  nation  membranes  in  polymer  electrolyte  fuel  cell,  J.  Elec¬ 
trochem.  Soc.  148  (3)  (2001)  A183-A188. 

[13]  F.  Barbir,  J.  Braun,  J.  Neutzler,  Effect  of  collector  plate  resistance 
on  fuel  cell  stack  performance,  in:  S.  Gottesfeld,  T.F.  Fuller  (Eds.), 
Proton  Conduction  Membrane  Fuel  Cells  II,  vol.  98-27,  The  Elec¬ 
trochemical  Society,  Pennington,  NJ,  1999,  pp.  400-406. 

[14]  B.  Andreaus,  A.J.  McEvoy,  G.G.  Scherer,  Analysis  of  performance 
losses  in  polymer  electrolyte  fuel  cells  at  high  current  densities 
by  impedance  spectroscopy,  Electrochim.  Acta  47  (2002)  2223- 
2229. 

[15]  F.  Barbir,  H.  Gorgun,  X.  Wang,  Relationship  between  pressure  drop 
and  cell  resistance  as  a  diagnostic  tool  for  PEM  fuel  cells,  J.  Power 
Sources  141  (1)  (2005)  96-101. 

[16]  F.  Barbir,  A.  Husar,  V.  Venkataraman,  Pressure  drop  as  a  diagnos¬ 
tic  tool  for  PEM  fuel  cells,  in:  The  2001  Electrochemical  Society 
Meeting,  San  Francisco,  CA,  2001. 

[17]  W.  He,  G.  Lin,  T.V.  Nguyen,  Diagnostic  tool  to  detect  electrode 
flooding  in  proton  exchange  membrane  fuel  cells,  AIChE  J.  49  (12) 
(2003)  3221-3228. 


